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An electrochemical STM is a rather innovative experimental method for measuring electron trans- 
port in single molecule nano-junctions, where for an atomistic interpretation of such experiments 
insight from density functional theory based calculations is needed. In this article we focus on co- 
herent electron transport as one of the possible transport mechanisms, where the explicit influence 
of the solvent on the transmission function has only rarely been included for neutral molecules in the 
literature and to our best knowledge never before for redox-active systems. Naively one would expect 
a rigid shift of the transmission function with the oxidation state with respect to the Fermi energy 
of the leads but we find more complex behaviour because of equilibrium charge transfer between the 
molecule and the leads in dependence on the number of counter ions, which we investigate in the 
context of electronegativity theory. For adjusting the oxidation state of the redox-active molecule we 
correct for self interaction effects by fixing the charge on the counterions, which in our calculations 
mimic the effect of the gate in the STM setup, with two competing methods, namely the generalized 
A SCF technique and screening with solvation shells. 4A comparison of the two approaches allows 
for a physical interpretation of the oxidation state dependence of electron transport in this junction. 



I. INTRODUCTION 

Most studies in the vibrant field of single-molecule elec- 
tronics focus on the low bias current flow through rather 
small benchmark molecules anchored to metal leads in ul- 
trahigh vacuum (UHV) at very low temperatures. Under 
those restrictions the underlying electron transport prob- 
lem is nowadays straightforwardly accessible to a compu- 
tational treatment with a nonequilibrium Green's func- 
tion (NEGF) approach^ in combination with a density 
functional theory (DFT) based description of the elec- 
tronic structure of the separate and combined compo- 
nents of the junction, namely the leads and the scatter- 
ing region 2 - - — . This method allows for an atomistic inter- 
pretation of associated UHV experiments on such bench- 
mark systems in a mechanical break-junction or scanning 
tunneling microscope (STM) setup^A thereby contribut- 
ing to a fundamental understanding of the dependence of 
the electronic conductance of the junction on the details 
of its structure within the boundary conditions of a low 
pressure and low temperature regime. 

For single-molecule junctions to be useful as molecular 
devices, however, their operability at room temperature 
is required and the presence of a solvent allows for elec- 
trochemical gating, which makes it possible to avoid the 
potentially destructive effect of the rather high local elec- 
tric fields, which otherwise would be needed for inducing 
a larger current^. Experimentally, these ambient condi- 
tions can be achieved with an electrochemical STMiSr— , 
where the nano junction is an integral part of an elec- 
trochemical cell and the investigated molecules usually 
have a redox-active center with an oxidation state which 
can be regulated via gating^. Depending on the setup 
as well as structural details of the system, two compet- 
ing electron transport mechanisms have to be considered 
for a theoretical description of such experiments, namely 
electron hopping which is a thermally induced two-step 



process and coherent tunneling which is the standard 
one-step phenomenon known from benchmark molecules 
without a redox-active center and relatively strongly cou- 
pled to metallic electrodes at temperatures close to 
K. In both cases an atomistic description of the process 
under electrochemical conditions provides a formidable 
challenge for a DFT based theory. While for the for- 
mer, the difficulty lies in a simplified and compact but 
nevertheless sufficiently accurate description of the nu- 
clear vibrations of the molecule and solvent which drive 
the electron flow, for the latter it becomes necessary to 
adjust the oxidation state of the redox active center in 
the scattering region and therefore deal with the issue of 
charge localization in a multi-component system. 

A correct description of localized charges is notoriously 
hard to achieve within a DFT framework, because the 
Coulomb and exchange parts of the interaction of an elec- 
tron with itself do not cancel out exactly in a standard 
Kohn-Sham (KS) Hamiltonian and the corresponding self 
interaction errors (SIE) result in an artificial tendency 
towards dclocalization^r— . As has been shown recently, 
both for a continuum solvation modcl^ and for an ex- 
plicit description of a periodic cell with its vacuum part 
filled up with H2O molecules^, a polar solvent has a 
screening effect on the Coulomb potential which reduces 
SIE and stabilizes localized charges within DFT. Another 
way to enforce localization is based on the generalized 
A SCF techniqu o 20 ' 21 , where an arbitrary integer value 
between and 2 for the occupation number of a partic- 
ular crystal eigenstate or linear combination of crystal 
orbitals can be defined as a boundary condition to the 
self-consistency cycles determining the electronic struc- 
ture of a given system. 

In our article we pursue both avenues for a 
study of the coherent electron transport through the 
Ru(PPli2)4(C2H4)2 bis(pyridylacctylydc) complex in 
Figure [TJ which we will often refer to as just "the Ru- 



2 




FIG. 1: Geometry of the Ru(PPIi2) 4 (C 2 H4)2 
bis(pyridylacetylyde) complex studied throughout our 
article bonded to ad-atoms on Au fee (111) surfaces 
within an aqueous solvent containing CI counterions. 



complex" in the following since it is the only system we 
investigate here and where for experiments in an aequous 
solution with chlorine counterions the oxidation state of 
the redox active ruthenium atom can be switched be- 
tween +11 and +III by varying the electrochemical poten- 
tial of the cell. We chose this particular system because it 
was used in previous conductance measurement :] 22 ' 23 as a 
monomer of chains - albeit with different anchor groups - 
where it was found that depending on the chain length ei- 
ther coherent transport or electron hopping is observed^ 2 -. 
In addition spectroscopic and quantum chemical studies 
on similar Ru complexes 2 ^ - — suggest that this molecular 
species offers the possibility to easily link two carbon- 
rich chains to each other for the formation of reversible 
redox systems 2 ^— with an optical transition highly de- 
pendent on the nature of the chai n 32 ' 33 , thereby serving 
as a starting point for the investigation of chains with 
multiple redox active centers 2 ^. In contrast to Refi^ we 
use pyridil groups as anchors to the leads because they 
provide peaks in the transmission function, which are 
narrow enough to assume that the oxidation state of the 
complex has an impact on the conductance but broad 
enough to avoid the Coulomb blockade regime 3 ^—. 

Although conductance calculations on redox- active 
complexes have been published befor o 37 ' 38 , we believe 
our article to be the first DFT based study of coher- 



ent electron transport through such a molecular complex 
which explicitly considers the aequous solvent and also 
investigates the influence of the oxidation state on the 
resulting transmission function. There have been pre- 
vious studies on the impact the solvent has on smaller 
benchmark molecules without a redox center— ~—, where 
some of the m 40 ' 42 have found a "chemical gating" effect, 
i.e. a shift in the transmission function induced by the 
surrounding molecules, which was explained by dipolc 
fields. We do not consider configurational fluctuations of 
the solvent molecules in our article, not only because of 
the high computational demands this would generate for 
our rather large junction but also because it would lead 
to fluctuations in the oxidation state of the Ru complex 
where a main aim in this work is to keep it fixed. 

The paper is organized as follows: In the next sec- 
tion we present transmission functions for the Ru com- 
plex at oxidation states of 0,+I,+II and +III where for 
mimicking the gate potential the counter charge is local- 
ized on a varying number of chlorine ions and the two 
ways of reducing SIE mentioned above, namely employ- 
ing the generalized A SCF technique and introducing 
H2O molecules explicitly as a solvent are used for this 
purpose. In Section [II] we also compare the shift in pro- 
jected molecular eigenvalues in dependence of the oxida- 
tion state which both methods produce. Section IIIII is 
devoted to the distribution of partial charges throughout 
the junction, which is a multi-component system in the 
sense that varying the oxidation state does not only af- 
fect the charge on the Ru-complex and counterions but 
also the gold leads and aqueous solvent can and do lose 
or gain fractions of electrons. For the analysis of this 
complex behaviour we start from cluster models within 
the simplified picture of electronegativity (EN) theory^ 3 -, 
and from their direct comparison with our full calcula- 
tions on the junctions represented by Figure [1] we derive 
the nature of the driving forces, which define the charge 
density distributions we observe. We conclude with a 
brief summary of our results. 



II. ELECTRON TRANSPORT CALCULATIONS 

All calculations of transmission probabilities T(E) in 
this article were performed within a NEGF-DFT frame- 
work^ - — with the GPAW cod e 44 ' 45 , where the core elec- 
trons are described with the projector augmented wave 
(PAW) method and the basis set for the KS wavefunc- 
tions can be optionally chosen to be either a real space 
grid or a linear combination of atomic orbitals (LCAO), 
and we opted for the latter on a double zeta level with 
polarisation functions (DZP) for all of our electron trans- 
port and most of our electronic structure calculations. 
The sampling of the potential energy term in the Hamil- 
tonian is always done on a real space grid when us- 
ing GPAW, where we chose 0.18 A for its spacing and 
a Perdew-Burke-Ernzerhof (PBE) 4 ^ parametrisation for 
the exchange-correlation (XC) functional throughout this 
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article. 

In this section we present our results for the trans- 
mission functions, where we highlight the methods we 
use to adjust the oxidation state of the Ru-complex and 
to cut down on the high computational demands of the 
NEGF-DFT procedure for this rather large system. In 
order to form a more direct connection between the elec- 
tron transport through the junction and the molecular 
orbitals (MOs) of the molecule, we also show their pro- 
jected energetic positions with respect to the Fermi level 
(E F ) of the leads by decoupling the basis functions situ- 
ated on the molecule from those on the Au surface and 
solvent in the transport Hamiltonian and perform a sub- 
sequent sub-diagonalization of its molecular par t 47 i 48 . 

A. Transmission functions 

Within NEGF the transmission function T(E) is de- 
fined by 

T(E) = Tr{G d T L G\Y R ) (1) 

where 

G d = (£ - tf d - £ L - E^r 1 (2) 

represents the Greens function of the device containing 
the self energy matrices ^l/r due to the left/right lead, 
T L / R = i(T, L / R — T^ L ^ R ) and the Hamiltonian matrix 
for the device region, which contains not only the Ru- 
complex but also 3-4 layers of the aligned Au surface on 
each side. Due to the rather large size of the central 
molecule (Figure [T]), we had to use gold slabs with a 
6x6 unit cell in the surface plane in order to ensure that 
neighbouring molecules do not interact. With the two 
Au ad-atoms directly coupling to the molecule (Figure 
[T]), the device region contains a total of 254 Au atoms 
in addition to the atoms of the complex itself and up to 
64 H2O molecules. As a consequence reached a size 
which was beyond our computational capabilities to be 
handled efficiently for electron transport calculations and 
therefore needed to be reduced as we describe in detail 
in the following paragraph. 

Since it is known that the solvent does not contribute 
to the peak structure in T(E)i£, but instead adds a 
base-line conductance with a rather small energy depen- 
dence^, we cut out the lines and rows indexing H 2 basis 
functions in the matrix H^, which we initially obtained 
from an electronic structure calculation for the full de- 
vice region. In a second effort towards memory reduction 
we cut out very high- and very low-lying MOs from 
after sub-diagonalizing it with respect to molecular basis 
functions, where we assumed that molecular eigenstates 
which arc further than 5 eV apart from Ej? would have 
no effect on the conductance or on the transmission func- 
tions on the much smaller energy range on which we show 
them. 



In the following wc discuss how to adjust the oxidation 
state of the Ru-complex in order to assess its influence on 
the transmission function and conductance of the junc- 
tion. For ensuring overall charge neutrality in the unit 
cell of our device region which is a necessity when apply- 
ing periodic boundary conditions for electronic structure 
calculations, the counter charge to the positively charged 
Ru-complex has to be an explicit part of the cell and in 
our system is represented by CI countcrions. There are 
two methods we exploit in this article to overcome the SI 
problem, which leads to an artificial derealization of oth- 
erwise localized charges in DFT: i) We make explicit use 
of the findings of other group o 18 ' 19 that a polar solvent, 
H2O in our case, stabilizes localized charges, because the 
solvation enthalpy and therefore also the total energy of 
the system become the more negative, the more point-like 
the charges on the solutes arc distributed; ii) and we also 
employ the generalized A SCF techniqu o 20 ' 21 ' 45 which 
has been previously used as a feature of GPAW for a cor- 
rect description of excitation processes in molecules ad- 
sorbed on surface s 19 ' 20 and of electron hopping between 
layers of oxide a 49 ' 50 . 

In practical terms the first scheme starts with the 
relaxation of the nuclear positions of the isolated Ru- 
complex towards the convergence criterium of 0.05 cV/A 
for the average force. Then we add one counterion with 
a fixed Ru-Cl distance of 7 A and embed the resulting 
system in a solvent shell of 46 molecules by making use 
of the graphical interface of the ghemical code^ 1 -, which 
places H2O molecules in the cell with a high degree of 
artificial translation symmetry. In a next step we relax 
the nuclei of the system now comprising the complex, 
the counterion and the solvation shell in order to create 
a more natural distribution of water molecules, where 
hydrogen bonds create a network structure, but we keep 
the Ru-Cl distance constant as a boundary condition for 
avoiding hybridization between the Ru-complex and the 
chlorine ion, which is statistically unlikely in nature but 
might happen in our relatively small unit cell. During 
this relaxation process we regularly probe the charge dis- 
tribution in the system. Once we achieve a one-electron 
charge on the Ru-complex, i.e. an oxidation state +1, wc 
stop the relaxation and align the whole system between 
two gold (111) surfaces with ad-atoms and the nitrogen 
of the pyridil anchors at a distance of 2.12 A for estab- 
lishing the direct electronic contact^. For this system 
we then calculate the transmission function as described 
above. 

In an attempt to adjust an oxidation state of +11 on 
the Ru-complex we add another CI atom again with a 
distance of 7 A to the Ru-position and now enlarge the 
solvation shell to a size of 58 water molecules. This 
enlargement of the size of our system which is at this 
stage not yet aligned to the Au-slab requires a repetion 
of the relaxation procedure described above but this time 
our relaxation converges below the given maximum force 
threshold without the achievement of a +11 charge on 
the complex. This behaviour can be explained by the 
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FIG. 2: Transmission function of the Ru-complex in 
different oxidation states adjusted with two different 
methods, i.e. a) solvent screening and b) ASCF. In 
both cases CI atoms were used as counterions, to 
extract the electrons from the Ru-complex, while 
ensuring charge neutrality within the cell. The kpoint 
sampling was performed on a 4x4x1 mesh in all cases. 



increase of the electronegativity of the Ru-complex with 
its oxidation state as we will do in some detail further 
below in Section MI Bl At this point we just state that 
we proceeded by adding another CI atom which enabled 
us to adjust +11 on the complex. Reaching an oxidation 
state of +III has proven to be impossible by applying the 
described method. We managed, however, to achieve a 
total charge of -III on five CI counterions with a solvent 
size of 63 water molecules, where a high percentage of 
the electrons are extracted from the solvent molecules, 
as will be further discussed in Section [ill Al Any further 
increase in the positive charge on the Ru-complex would 
have required a much larger solvation shell, which was 
ruled out by the corresponding computational costs also 
including a substantial augmentation of the gold slab. 

In our second approach based on the generalized A 
SCF method, we make use of its flexibility to define 
the spatial expansion of an orbital enforced to contain 
an electron as an arbitrary linear combination of Bloch 



TABLE I: Conductance of the Ru-complex in the setup 
of Figure [1] in different oxidation states as calculated by 

NEGF-DFT. The values are given in conductance 
quanta Go, where the value of the complex with 

oxidation state was calculated to be 1.6T0 -5 Gn. 



state a 20 ' 21 . By extracting one electron from the system 
and inserting it into a predefined orbital in the begin- 
ning of every iteration step, the self consistency cycle 
progresses as usual but with the electron density of this 
particular orbital as a contribution to the external poten- 
tial. The same process can be applied for two or three 
electrons. In this way we can fix the electron occupation 
on the counterions manually, which solves the self inter- 
action problem implicitly and makes this method ideal 
for charge localization as needed in the present work. 
When applying this technique we chose the nuclear posi- 
tions relaxed for the neutral complex aligned betwen the 
gold surfaces, where the respective number of counterions 
were added, each with one supplementary electron con- 
strained to completely fill its p orbital. This procedure 
also had the benign consequence that the calculation of 
T(E) was reduced significantly in terms of computational 
demand, because we do not need an explicit solvent here 
and therefore do not have to remove the respective states 
from the transport Hamiltonian. 

Figure [2] shows the transmission function calculated 
for the Ru-complex in different oxidation states adjusted 
with the two methods described above. For a comparison 
with experiments only the oxidation states +11 and +III 
are relevan t 22 i 23 , because only they correspond to stable 
compounds, but we nevertheless cover and +1 in our 
article as well where we aim at a systematic investigation 
of the dependence of electron transport properties on the 
charge distribution in the system. Naively one would ex- 
pect that the +III charged complex has a higher conduc- 
tance than the +11 charged one due to its half filled MO 
at the Fermi Level. As shown in Table U this is actually 
the case when the charge is localized by ASCF but not 
when we use the method where we stabilize charges with 
solvation shells. Part of the reason for this disagreement 
can be explained by looking at Figure [5^,, where we find 
that the incompleteness in decoupling the H2O orbitals 
from the transport Hamiltonian -where one has to con- 
cede that LCAO basis functions located on specific atoms 
also contribute to the description of its surrounding- cre- 
ates a "transmission baseline" which fits to the behaviour 
previously investigated in theoretical studies of the con- 
ductance of water—. This contamination issue is also the 
reason for the artificial peak, which can be seen slightly 
above the Fermi energy for the 5C1/63 water junction in 
Figure [5^. 
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This water baseline for the electron transmission is ob- 
viously absent in the ASCF calculations where we find a 
continuous rise in conductance with an increase in oxida- 
tion state (Table [TJ coming from the steady progression 
of a peak towards ~Ep from below (Figure [2J)) . As we 
will point out in the next subsection, where we focus 
on the molecular eigenstates, this double-peak is related 
to the HOMO and HOMO-1 of the Ru-complex which 
are similar in their symmetry and spatial resolution. It 
is the common understanding of such systems, however, 
that the HOMO should cross the Fermi level when mov- 
ing from +1 to +11 and the half-filled orbital at +III 
should be the HOMO-1 which is clearly not the case for 
the transmission functions calculated with either of our 
two techniques and displayed in Figure [2] For a further 
analysis of our findings a detailed study of the charge 
density distribution in the junction is needed which we 
will provide in Section IIII Al 



B. Projection of MO eigenenergies 

In order to understand the peak structure in Figure 
[5] in more detail we now study the electronic structure 
of the junction by investigating the electronic states of 
the device in terms of the molecular eigenenergies and 
their shape. Since the coupling of the Ru-complex to the 
Gold (111) surface leads to a hybridisation of the respec- 
tive electronic states, it is necessary for the extraction 
of the molecular eigenvalues of the Hamiltonian matrix 
to eliminate their coupling to the surface states so that 
the MOs on the Ru-complex become isolated. In a sub- 
diagonalization procedure, this is done by setting all the 
matrix elements of the transport Hamiltonian to zero, 
except the ones which correspond to the basis states for 
the molecule. 

The MO-eigenvalue distributions obtained in this way 
are shown in Figures [3] While the left panel displays the 
eigenvalues of the system prior to its attachment to the 
gold surfaces, the right panel shows the eigenvalues after 
dehybridisation of a Hamiltonian which was based on the 
full electronic density of the device region including its 
Au part. In Figure [3J the MOs on the counterions and 
the solvent are plotted explicitly to investigate the influ- 
ence of the charge transfer in the whole system. Naively 
one would expect the following three-step behaviour of 
the MO eigenvalues of the Ru-complex without contact 
to the Gold surface: 1) Charging the molecule to its +1 
state extracts one electron from the complex's HOMO, 
which leads to a SOMO, that by definition is situated at 
the Fermi level. 2) When in a next step another electron 
is extracted from the molecule the former HOMO is emp- 
tied, which must lead to an energy above the Fermi level. 
3) At an oxidation state of +III the former HOMO-1 is 
now half filled and defining the Fermi level. In this ideal 
picture the CI p-orbitals are always fully filled, with their 
HOMO lying below the Fermi level. As can be seen in 
the left panel of Figure [3] these expectations are well met 



when the charge on the countcrion is defined by ASCF. 
This also applies for the +1 and +11 charged cases, when 
the oxidation state was defined by counter ions and sol- 
vent, while here only the +III system is not characterized 
by a half filled HOMO-1. 

The chlorine p-states are always situated below the 
Fermi energy in the left panel of Figure which suggests 
that they are indeed fully filled, but at higher oxidiation 
states they shift closer to Ep. Since more counterions 
were used then electrons had to be extracted from the 
complex the p-states are then only partially filled by def- 
inition for oxidation states +11 and +III. It has to be 
noted that also water takes part in this multi compo- 
nent system, whose role as a charge donor or acceptor 
should not be disregarded. Without the gold slab the 
water states define and therefore some of them have 
to be partially filled orbitals. 

The electronic coupling of the complex to the Au sur- 
face changes the energy landscape of the system in several 
ways as can be seen in the right panel of Figure [3j Due 
to its metallic character and the large number of gold 
atoms, the Fermi level is now defined mostly by the gold 
states. As a result already the eigenvalues of the com- 
plex without counterions i.e. the uncharged complex, 
are shifted to higher energies relative to their new refer- 
ence point. Raising the oxidation state of the complex 
now leads to an eigenvalue alignment^— which is signif- 
icantly different from the isolated complex. The HOMO 
of the Ru-complex now always lies distinctively below the 
Fermi energy. While gold creates a new reference energy 
for the molecular eigenstates, it also plays the role of an 
electron donor or acceptor, meaning, that it can accept 
charge from both the complex and the counterion/solvent 
system. We also note in this context, that for pyridil an- 
chors on Au Pauli repulsion leads to an electron depletion 
on the comple»^. As a consequence the HOMO of the 
Ru-complex, does not cross the Fermi energy for any ox- 
idation state in the right panel of Figure [3J i.e. it never 
becomes empty entirely. Assuming that the counterions 
are the only charge extracting part in the system and 
taking into account that their p-states seem to be fully 
occupied judged by their distance from the Fermi level in 
the right panel of Figure [3j the positive charge which is 
not localized in the complex has to be found in the gold 
surface or the water molecules, where when comparing 
the two panels of Figure |3] a significant change in the 
cigenencrgy spectrum of the water molecules due to the 
presence of the gold surface can be observed. Since we 
are dealing with a four component system for the charge 
transfer here (Ru-complcx/gold-surface/chlorine/water) 
a detailed charge density distribution analysis has to be 
performed for its interpretation, which will be the topic 
of the following section. 

By visual inspection of the shape of the relevant or- 
bitals for coherent transport, i.e. the HOMO and the 
HOMO-1 of the Ru-complex, which we show as insets 
in Figure O we find that both MO's are characterized 
by a conjugated 7r-system, which is delocalized over the 
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FIG. 3: MO eigenvalue spectrum of the device region without (left) and with (right) coupling to the gold slab, where 
for (b) the MOs are calculated by decoupling the basis functions localized on the molecule from that of the surface 
states with a subdiagonalization of the transport Hamiltonian. Each panel is divided into several sections, namely 
the electronic eigenvalues of the uncharged Ru-complex, followed by its MOs with the charge constrained by ASCF 
and its MOs from the solvent screening method (the spatial shape of the HOMO and HOMO-1 are given as insets 
for both methods), where also the eigenvalues of the chlorine ions (green) and water molecules (blue) are shown. 



whole bridge of the complex and their energies explain 
the double-peak structure in the transmission function 
in Figure [3J There are, however, some differences be- 
tween the two MOs. While the HOMO in the left panel 
of Figure [3] has a high localization at the interface re- 
gion, the HOMO-1 does not, which should diminish its 
coupling strcngh to a possible surface in the latter case. 
Comparing these MO's in Figure [3] shows that the actual 
contact of the Ru-complex to the Au slabs changes their 
energetical hierarchy with the HOMO-1 changing its role 
to being the HOMO now and resulting in a narrower peak 
in the transmission functions in Figure [2] than the other 
MO due to its weaker coupling to the surface. 



III. EQUILIBRIUM CHARGE DISTRIBUTION 

The key for understanding the peak positions in the 
transmission function and the Fermi level alignment of 
the corresponding MOs lies in understanding zero bias 
charge transfer as has been shown in Refs^r— . The 
present case, however, is more difficult because here we 
have to deal with a four component system (Ru-complex, 
CI, solvent, leads) and with self interaction effects. In 
the following we argue that the SI problem is overcome 
by both methods, namely the screening of the point 
charge on the counterion with a solvent and the con- 
straint of orbital occupation with ASCF, which allows 
us to reduce the complexity of the investigated system 
by replacing the chlorine ions and solvent by an external 



charge, within the context of electronegativity theory for 
our analysis. Throughout this section we use the Bader 
metho d 56 ' 57 for the definition of the electronic charges 
belonging to particular nuclei, which attributes fractions 
of the electron density by searching for minima on a grid. 
This division of the cell in small segments makes it possi- 
ble to define so called Bader volumina, and by integration 
over them the charge on a certain atom can be defined. 



A. Charge distribution in multi-component 

systems 

Tabic |TT] shows the charge on the Ru-complex in depen- 
dence on the basis set and the number of k-points, for the 
calculations with and without a Au slab in the cell. For 
all cases the geometries are optimized to reproduce the 
oxidations states +I,+II and +III respectively, so that 
the chlorine atoms and solvent molecules contain the re- 
spective counter charge as defined by the Bader analy- 
sis. The charge values on the Ru-complex without the 
Au surface correspond to the requested oxidation state 
rather precisely, and for the +1 and the +11 state both 
applied charge localization methods show no significant 
difference in their performance and basis set dependence, 
while the +111 state can only be reached with the ASCF 
occupation constraint. An explanation for the limits of 
the localization method with the solvent screening lies 
in the charge distribution between the CI atoms and the 
water. For the first two oxidation states the solvent plays 
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FIG. 4: Charge density difference between the coupled 
system and the isolated complex and gold slab, where 
pseudo densities in terms of the PAW formalism have 
been used for the densities, in order to eliminate 
artificial peaks near the nuclei. 



the role of an electron acceptor, meaning that the elec- 
trons which are extracted from the complex are situated 
at the CI atoms and the surrounding water molecules. 
Summing up the charge on the water molecules therefore 
leads to a negative charge on the solvent i.e. an excess of 
electrons, which originate from a corresponding electron 
loss on the Ru-complex. With five counterions in the 
solvent we find a contrary behaviour for the +111 case, 
i.e. while the CI atoms still exhibit an electron surplus of 
2.8 electrons with a LCAO basis and even 2.95 electrons 
with a grid basis, the solvent molecules now play the role 
of an electron donor with a loss of 1.10 electrons with a 
LCAO and 1.05 electrons with a grid basis. 

Adding the metal surface to the simulation cell changes 
the charge distribution drastically. While in the analysis 
above only the Ru-complex lost electrons to the counte- 
rion/solvent system, now also the metal surface can con- 
tribute in the same way. As can be seen in Table [TTTI even 
the complex with oxidation state loses electrons to the 
Au surface which we attribute to Pauli repulsion^. In 
agreement with that assessment, Figure |4] shows that the 
charge transfer happens mostly at the interface, while the 
rest of the junction is not contributing to it in a signifi- 
cant way. This behaviour poses an accuracy challenge for 
the Bader analysis, which is not best suited to describe 
diffuse states and therefore of the numbers in table HT1 
and IIII1 have an uncertainity in the range of 0.01-0.05 
electrons. 

The charge rearrangment at the interface between sur- 
face and molecule due to Pauli repulsion is independent 
of the chosen method to localize the electrons, but basis 
set dependent. This is particularly the case when ASCF 
is applied, which due to its nature of defining the selected 
orbital in terms of the basis functions of the system is by 
way of construction strongly dependent on the quality of 
this basis itself. The results from the solvent screening 
method on the other hand have proven to be rather sta- 
ble in terms of the used basis. By the introduction of 
a kpoint sampling the lack of accuracy in the definition 
of the constrained state in ASCF can either cancel out 
by the related averaging or sum up, which can be seen in 
Table ILTl where we compare the charge on the complex for 



a finite number of kpoints with that from a single k-point 
calculation for the ASCF approach. 

For both methods for charge localization the gold bulk 
tends to absorb increasing amounts of the positive charge, 
as shown in table IIII1 which explains why the HOMO 
peaks in the transmission functions shown in Figure [2] do 
not cross the Fermi Level for higher oxidation states in 
spite of a counter charge higher than 2 electrons. It is 
a delicate question wether this is due to SI artefacts in 
the calulations or a realistic result for the investigated 
system. While we deal explicitly with the SI error for 
the charge localisation on atoms, the charge distribution 
between the Ru-complex and the gold slab is not nec- 
essarily strongly localised anywhere, since the Ru atom 
is embedded into the complex by rather strong covalent 
bonds with its ligands. Also the electronic coupling at the 
interface is of intermediate strength, as indicated by the 
rather broad peak shape in the transmission functions. 
This does not contradict with the fact, that the bonding 
between the pyridil anchor group and gold atom is rather 
weak^, because in the case of Pauli repulsion the cou- 
pling with filled MO's produces bonding and antibonding 
states^. So the charge distribution we find in table IIIII 
could be physically correct, although it is not what one 
would attribute to the system when writing down its re- 
dox equations. For investigating the issue whether the 
charge distribution between the gold slab and the Ru- 
complex is realistic further, we employ electronegativity 
theory in the next subsection. 



B. Electronegativity theory 

In order to enable a view on our results from a different 
perspective we now analyze the junction in terms of elec- 
tronegativity theory following the concepts of Parr and 
Pearson^. The key quantities in this approach are the 
electronegativity \i and the hardness v 1 where the first is 
based on Mulliken's definition of electronegativity^, i.e. 



- f—\ - I + A 



(3) 



and the latter is defined as 



1 fd 2 E\ 

2 I ON 2 I 



I- A 



(4) 



with I being the ionisation potential, calculated as the 
total energy difference of the N and N-l system, and A 
the electron affinity, defined as E(N + 1) - E(N). 

When two different systems are brought into contact 
the charge transfer from one to the other can be calcu- 
lated as 



AN 



^2 - Mi 
2(^i + 1*2) 



(5) 
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ASCF 


Solvent 




1C1 


2C1 


3C1 


1C1/46 H 2 


3C1/58 H 2 


5C1/63 H 2 


Without Gold surface 


LCAO 1 kpoint 


-0.97 


-1.94 


-2.89 


-0.98 


-2.04 


-1.68 


Grid basis (h=0.15) 


-0.96 


-1.99 


-2.81 


-0.99 


-2.12 


-1.90 


Grid basis (h=0.18) 


-0.94 


-1.88 


-2.84 


-0.99 


-2.12 


-1.90 


Complex aligned to Gold bulk 


LCAO 8 kpoints 


-0.68 


-1.19 


-1.70 


-0.90 


-1.26 


-1.41 


LCAO 1 kpoint 


-0.73 


-1.33 


-1.97 


-0.89 


-1.24 


-1.41 


Grid basis (h=0.15) 


-0.72 


-1.42 


-2.26 


-0.84 


-1.22 


-1.37 


Grid basis (h=0.18) 


-0.85 


-1.51 


-2.23 


-0.88 


-1.22 


-1.43 



TABLE II: Basis set and k-point sampling dependence of the electron loss on the Ru-complex, for different 
quantities of counterfoils and solvent molecules. The two scenarios with and without the gold surface as part of the 

system are both displayed. 







Au 


H 2 


CI 


H 2 0+C1 


Ru complex 


Ru 




Neutral complex 


0.39 








-0.43 


-0.21 




1C1 


-0.16 




0.89 


0.89 


-0.725 


-0.25 


ASCF 


2C1 


-0.49 




1.80 


1.80 


-1.34 


-0.285 




3C1 


-0.69 




2.63 


2.63 


-1.97 


-0.29 




1C1/46 H 2 


-0.21 


0.37 


0.71 


1.08 


-0.90 


-0.24 


Solvent 


3C1/58 H 2 


-0.37 


-0.42 


1.95 


1.53 


-1.24 


-0.27 




5C1/64 H 2 


-0.58 


-0.69 


2.64 


1.95 


-1.40 


-0.37 



TABLE III: Distribution of the partial charges in the system including the Au slab as calculated with the ASCF 

and counterion/solvent schemes. 



where both the electronegativities and hardnesses of 
the separate components contribute to the regulation of 
the charge transfer—. 

The ionisation potential I and the electron affinity A 
are commonly defined for the neutral state of the indi- 
vidual subsystems, but as shown by Balbas et ali£°-, their 
role of defining the electronegativity and hardness is also 
valid for ions, which allows us to describe also the charge 
distribution in the junction at higher oxidation states in 
terms of EN theory. 

As discussed in the previous sections we fixed the 
charges on the counterions and solvent manually, and 
therefore we are mostly interested in understanding the 
charge distribution between gold and the Ru-complex, 
within the framework of electronegativity theory. For 
this purpose we adjust the oxidation states by defining 
an external charge q for the calculations on the subsys- 
tems, as indicated in Egns [3] and [4] as index of the partial 
derivatives, which is unusual, but not in contradiction to 
the basic assumptions of EN theory. It requires, however, 
some statistics for taking into account all possible start- 
ing points for the charge transfer. In the case without ex- 
ternal charge only one such initial electron configuration 
of the components has to be dealt with, i.e. neutral gold 
slab and neutral complex. Raising the external charge 



to +l\e\ already allows for two different starting points 
for the charge transfer, namely Ru-complex +1 /Au° and 
Ru-complex /Au +1 . In principle the calculation of AN 
for both should lead to identical predictions for the fi- 
nal charge distribution in the composite system with 
an oxidation state adjusted to +1, but imperfections of 
our DFT based total energy calculations such as SI er- 
rors and the approximative nature of the XC functional 
lead to deviations, as shown in table IIVI Averaging AN 
over all possible integer configurations should provide an 
improvement with regard to such errors. Hereby spe- 
cial emphasis has to be put on the reference point for 
AN, i.e. the subsystem with the index 1 in equation [5] 
Since the Ru-complex and the gold slab enter this equa- 
tion at charged states, the calculated fii,Vi and AN are 
also referring to these charged states. In order to ob- 
tain the change of electrons relative to the neutral sub- 
systems therefore the related integer charges have to be 
subtracted. 

For understanding the role of the size of the gold slab 
for the charge distribution we model the gold component 
in our EN theory analysis with clusters of different sizes, 
starting from the adatom and reaching up to the full gold 
surface used in the junction, as shown in Figure [SJ where 
we computed the electronegativities and hardnesses for 
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FIG. 5: Electron loss on the complex when an external 

charge up to +3|e| is applied. Panel (a) shows the 
values predicted from electronegativity theory (red) and 
from calculations where the complex is coupled to a 
gold cluster and the composite system was analyzed 
with the Bader analysis (black). In panel (b) the AN 
values from these two sets of model calculations are 
compared with periodic calculations (green line and 
blue crosses) of the device region, where the external 
charge was imposed as counter charge localized on CI 
ions and on solvent molecules. 



oxidation states from to +III for each cluster size in a 
setup without periodic boundary conditions and calcu- 
lated AN averaged over initial electron configurations as 
described above. In Table IIVI we show the related sta- 
tistical spread for the smallest and largest of our cluster 
sizes. Although we find that the deviations increase both 
with the order of the oxidation state and the size of the 
Au cluster, their overall values are reasonably small in- 
dicating that our predictions for AN from EN theory are 
not particularly limited in their accuracy by SIE or our 
choice of XC- functional. 

For building a bridge between the predictions for the 
charge distribution from EN theory and the actual ones 
we find in the periodic systems we use as device regions in 
the transport calculations, we also performed cluster cal- 
culations, containing both subsystems. The charge distri- 
bution in the resulting cluster cells were analyzed accord- 
ing to Bade r 56 ' 57 , where we imposed external charges for 



# of gold 
atoms 




254 




Starting charges 



Ru-complex 



+ 1 








Au-cluster 





+ 1 





+ 1 
+2 





+ 1 
+2 
+3 





+ 1 





+1 

+2 





+ 1 
+2 
+3 



AN 



-0.69 
-0.64 



-1.28 
-1.29 
-0.95 



-1.82 
-1.88 
-1.78 
-1.55 



-0.24 
-0.34 



-0.52 
-0.50 
-0.50 



-0.54 
-0.78 
-0.76 
-0.66 



AN 



-0.66 



-1.17 



-1.76 



-0.29 



-0.51 



-0.69 



TABLE IV: Illustration of the statistics in our EN 
theory predictions for charged states, which arises from 
the possibility of different initial charge configurations 
on the subsystems before they are brought into contact. 
The point of reference for AN in this table is the 
Ru-complex in its oxidation state 0. 



varying the oxidation state as we did for the subsystems 
for the EN predictions. 

The appeal of this intermediate step towards the peri- 
odic system calculation is that it allows us to distinguish 
between effects which come from electronegativity differ- 
ences of the components, others which find their origin 
in the spatial polarisation of the subsystem, when they 
are actually brought into contact in a given geometr y 54 ! 58 
and finally those related to the particular method we em- 
ploy for adjusting the oxidation state. 

As shown in Figure [SJt the results from the EN pre- 
diction and the Bader analysis of the composite systems 
in dependence on the Au cluster size differ. In this com- 
parison when we apply EN theory, the charge transfer is 
slightly underestimated for an external charge q=0 |e|. 
Raising q to finite values leads to an overestimation of 
AN with respect to the Bader analysis for the composite 
system. The deviation at high external charges is small 
for less then four gold atoms on both junction sides, but 
increases with the Au cluster size. 

Figure [SJj puts a different perspective on these qualita- 
tive differences when we compare the charge distributions 
obtained from both EN theory and the Bader analysis 
of the cluster calculations with the results for the peri- 
odic device region (Table H|. While in the latter case the 
charge on the molecule increases almost linearly with the 
counter charge, the two sets of model do not predict this 
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FIG. 6: Electron density difference between the 
Ru-complex in oxidation states +III and 0, where the 
oxidation was adjusted by an external charge (upper 
panel) or CI ions , with the negative counter charge 
localized by ASCF (lower panel). In both cases we show 

results from cluster calculations with an isovalue 
threshold of 2*10 -4 e, where a loss of electronic charge 
is depicted in blue and a gain in red colour. 



sign situated around the cluster. As a consequence a lo- 
cal coulomb attraction term makes due to its distance de- 
pendence a localization of the positive charge on the Ru- 
complex and the Au surface rather than the bulk regions 
more favourable. Figure [6] shows the charge density dif- 
ference between the highest charged state and the neutral 
junction, for the oxidation state defined by an external 
charge (upper panel) and by chlorine atoms with charge 
localization enforced by ASCF (lower panel). Without 
counterions the introduced positive charge is distributed 
almost completely over the gold atoms in the leads. Due 
to the non periodic setup of the cell the charges propa- 
gate to the outward pointing surfaces of the gold because 
of the mutual repulsion of the fractional positive charges. 
If the oxidation state is defined by chlorine counterions 
on the other hand, the introduced positive charge in the 
junction is mostly localized on the Ru-complex and the 
lead surface for the reasons just explained. Fractions of 
positive charge are, however, still localized on the outer 
parts of the gold bulk, because they are not hindered by 
the presence of a neighbour cell in a non periodic setup 
and Figure [SJd shows that therefore periodic boundary 
conditions even increase the positive charge on the Ru- 
complex region. On the basis of this analysis we be- 
lieve that the charge distributions presented in Table Hill 
which explain the qualitative behaviour of conductance 
and transmission in dependance on the oxidation state 
are physically reasonable. 



behaviour. The reason for this finding can be found in 
the detail of the oxidation state definition, where for the 
model calculations an external charge is imposed, which 
is distributed homogeneously over the subsystem. For 
254 gold atoms the introduced charge delocalizes over 
all the atoms in the cluster leading to just a minor in- 
crease of its electronegativity at higher charged states, 
which is a consequence of the hardness, as the derivative 
of the electronegativity (see equation [3] and [4]) becoming 
smaller with cluster size. On the other hand the energy 
needed to extract an electron from the much smaller Ru- 
complex increases strongly with its oxidation state com- 
pared to the gold. As a consequence the external charge 
is mostly absorbed by the Au cluster, leading to rather 
modest charging of the molecule with an increasing exter- 
nal charge. The effect is more pronounced in the calcula- 
tions for the composite junction with an external charge, 
where AN in Figure [S}) suggests, that the direct contact 
between the gold cluster and the Ru-complex enhances 
the hardness-driven charge localization on the gold slab 
when compared with EN theory predictions. 

If we adjust the oxidation state also in the composite 
cluster calculations in the same way we did for the peri- 
odic cells, however, a counter charge is localized on the 
chlorine and the situation changes, as can be seen from 
the blue curve in Figure [SJd. Instead of a globally defined 
external charge we now have point charges of opposite 



IV. SUMMARY 

The aim of this article was the description of coher- 
ent electron transport through a single molecule junc- 
tion containing a redox active center with an emphasis on 
its oxidation state, a scenario which to our best knowl- 
edge has never been studied on an ab initio level be- 
fore. A correct description of this oxidation state within 
DFT is essential and we applied two independent meth- 
ods for correcting the self interaction error, namely sol- 
vent screening and ASCF, where for both charge is lo- 
calized on counterions, which act similar to a gate in an 
electrochemical STM setup. We found that the actual 
charge on the Ru complex is smaller than its oxidation 
state, when it is coupled to a gold surface, which might 
indeed be realistic since some of the charge can be ab- 
sorbed by the leads. In order to investigate this issue 
we made predictions for model systems of varying size 
of the gold component within electronegativity theory, 
which we supplemented with cluster calculations. This 
analysis led us to the conclusion, that a large part of the 
charge should indeed be absorbed by the leads, but half 
of it remains on the complex due to coulomb attraction, 
where the vicinity of the localized charges on the counte- 
rions has a stabilizing effect. Our analysis led us to the 
conclusion that the charge distributions we find in our 
calculations for the device region are realistic in physical 
terms. 
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